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Abstract. Recent experimental results from the Trident laser confirm the im- 
portance of kinetic effects in determining laser reflectivities at high intensities. 
Examples observed include scattering from low frequency electron acoustic waves 
(EAWs), and the first few stages of a cascade towards turbulence through the 
Langmuir decay instability. Interpretive and predictive computational capability 
in this area is assisted by the development of Vlasov codes, which offer high ve- 
locity space resolution in high energy regions of particle phase space, and do not 
require analytical pre-processing of the fundamental equations. A direct Vlasov 
solver, capable of resolving these kinetic processes, is used here to address fun- 
damental aspects of the existence and stability of the electron acoustic wave, 
together with its collective scattering properties. These simulations are extended 
to realistic laser and plasma parameters characteristic of single hot-spot exper- 
iments. Results are in qualitative agreement with experiments displaying both 
stimulated Raman and stimulated electron acoustic scattering. The amplitude of 
simulated EAWs is greater than that observed experimentally, and is accompa- 
nied by a higher phase velocity. These minor differences can be attributed to the 
limitations of a one-dimensional collisionless model. 
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1. Introduction 

Recent single hot-spot experiments using the Trident laser facility ^1 E] identified 
backscatter which resembles stimulated Raman scattering (SRS) but, importantly, 
is due to stimulated scattering from an electron plasma mode whose frequency is 
significantly below the plasma frequency. This mode was identified as the electron 
acoustic wave (EAW), an undamped electron mode absent from fluid descriptions of 
the plasma. In order to survive and propagate, the EAW requires a non-Maxwellian 
electron velocity distribution, flattened at the phase velocity of the wave, without 
which it would be critically damped. We demonstrate that in laser-plasma conditions 
which prohibit SRS, an incident electromagnetic wave can drive an initially critically 
damped electron plasma wave to sufficient amplitude that it can trap electrons, which 
then support it and allow it to propagate undamped, forming an EAW. Alternatively, 
flattened distributions supporting an EAW may be created by the trapping of 
electrons during the saturation of laser-plasma instabilities, such as conventional 
SRS in the case where plasma conditions permit it. The physical and mathematical 
characteristics of the EAW are described here, and the mode is then simulated using an 
Eulerian Vlasov code to demonstrate its undamped nature. The implications for laser 
plasma interactions are examined, including the Langmuir decay instability (LDI) and 
stimulated electron acoustic scattering (SEAS). 

In particular, we model a scenario where the Langmuir wave excited by stimulated 
Raman scattering grows to an amplitude sufficient to cause a significant local flattening 
of the electron velocity distribution function. This local flattening is formally 
equivalent, as we shall see, to the creation of a small drifting beam population, 
and the EAW can be considered to be an electrostatic mode supported primarily 
by this beam population. Stimulated scattering of the incident laser light from the 
EAW then becomes possible (SEAS), by analogy with SRS. We show that numerical 
implementation of a model combining a kinetic (Vlasov) description of longitudinal 
dynamics with a fluid description of transverse dynamics successfully captures the key 
physics. 

In this paper we explore questions which, in addition to assisting the 
interpretation of present and possible future experiments in laser-plasma interactions, 
raise several interesting theoretical considerations. First, the literature on the theory 
and modelling of the EAW, which displays paradoxical qualities, is sparse, and for 
this reason we offer an ah initio treatment. Second, SEAS provides interesting points 
of contact between the phenomenology of energetic particle populations in laser- 
plasma interactions and in magnetically conflned fusion plasmas. These points of 
contact include: the existence of modes supported primarily by the energetic particle 
population; the independent role of such modes in coupling the plasma to external 
drivers; their role as a channel for energy transfer within the plasma, for example 
through a turbulent (or at least nonlinear) cascade; and their diagnostic potential. 
See, for example Refs. jHlElElE] and references therein. Third, there is the question 
of what level of theoretical description (kinetic, fluid, etc.) best captures the key 
physics while also enabling effective numerical simulation. 

2. The Electron Acoustic Wave 

The possible existence of plasma waves at frequencies significantly below the electron 
plasma frequency was first identified by Stix[7] , although it was expected that Landau 
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damping in this regime would prohibit their formation. Later work[5| 1^1 IIUI ITT] 
showed that EAWs can indeed exist, supported by a population of trapped electrons. 
Derivation of the plasma dispersion relations, based on a two-fluid treatment, yields 
high frequency Langmuir waves and low frequency ion acoustic waves. A more 
complete kinetic treatment, based on linearising the Vlasov equation, shows these 
modes to be damped. This is due to Landau damping, a purely kinetic effect which 
(in one dimension for simplicity) requires that dvf\vp < 0, where / is the particle 
distribution function and Vp the phase velocity of the wave. However if dvf\v^ ~ 
then the wave may be undamped and, as we shall see, this is the case for the EAW. 

In order to construct a dispersion relation for the EAW let us first consider a 
Maxwellian distribution with characteristic velocity vt that has a flattened region at 
V = Vp; we first take the limit where the width of the flattened region, in velocity 
space, tends to zero while dvf\vp = 0. The integral in the Landau dispersion relation 

e(w, fc) = 1 - — — 2 / — — -dv (1 



can be written 



,v-uj/{kvT) J V - uj/{kvT) 

By construction, the second term is zero, leaving only the principal value integral 
Evaluating this integral gives the dispersion relation for EAWs in the linear limit: 



fc^A^ + l-V2^Daw(-^).0 (3) 

Daw it) = exp (-t^) / exp(u )du (4) 
Jo 

where Eq.Q) is Dawson's integral ^21. Equation (O, which can be evaluated 
numerically, thus gives the dispersion relation for undamped plasma waves in the limit 
of vanishing amplitude. The dispersion curve, labelled Ai; = 0, is shown in Fig. ^ 
where there are two distinct branches. The upper branch corresponds to an undamped 
form of the Langmuir wave and the lower branch, with uj < ujpe, corresponds to the 
electron acoustic wave. 

Similar analysis can be carried out for distribution functions that have a finite 
flattened region of width Au proportional to the wave amplitude, defined by 

/ = /o + /i (5) 
where /q represents a Maxwellian distribution, and flattening a,t v = Vp is provided by 

h{v) = a./ok, {v - Vp) exp {^ '^""^J"^^ ) ■ (6) 

Here the width of the flattened region is related to the EAW amplitude by assuming 
that electrons trapped in the wave potential can be considered as simple harmonic 
oscillators. Equating the potential and kinetic energies of a trapped electron then gives 
the width of the structure in velocity space, which dictates Aw in Eq.@. Equivalently 
one may calculate the maximum amplitude of an EAW which can be supported by a 
given flattened velocity distribution. We flnd, given Eq.®, that Aw = y'e^/rTT^ where 
(j) is the wave potential. Assuming a sinusoidal waveform this gives Au — eEo/kme, 
where Eq is the wave amplitude; we note that kAv equates to the nonlinear bounce 




Figure 1. Dispersion relations calculated from Eq.(3J for undamped plasma 
modes in the linear limit (where the width Av of the flattened region, in units 
of thermal velocity, tends to zero) and for two nonlinear cases {Av = 0.5, 1.0 in 
Eq.jnj). Wavenumbers are normalised to and frequencies to uipe- The lower 
branch represents the electron acoustic wave, which for low k follows lu = 1.35fc. 
The upper branch represents an undamped form of the Langmuir mode. 



frequency of resonant particles, as defined by Eq.(27) of Ref . , for example. With 
/i defined in this way, Eq. ^ can be solved numerically to give a family of dispersion 
curves inside the ideal, infinitesimal amplitude case, as shown in Fig. ^ The dispersion 
relation of the EAW is thus linked to the wave amplitude. 

The EAW appears, at first, to be unphysical. In particular its low frequency 
is a characteristic not expected of electron plasma waves supported by distribution 
functions arbitrarily close to Maxwellian, in which ion dynamics play no role. Some 
physical understanding, for the case of small k, can be gained by considering a locally 
flattened distribution to be a superposition of a background Maxwellian population 
and a smaller drifting population, as shown in Fig[5J 

In the frame of reference of the background population, oscillations at the electron 
plasma frequency supported by the drifting population will be Doppler shifted such 

1 /2 

that u! = uj' + kv where lu' = (e^n2/eoWe) is the electron plasma frequency for this 
second population, whose density is n2- In the limit where n2 tends to zero, we recover 
u! ~ vk. This dispersion relation is linear in k for small k < 0.2X~jj^, and implies a 
frequency below the plasma frequency. While this interpretation does not give the 
exact value for the phase velocity of the EAW in the linear limit (vp w l.SbvTe, as in 
Fig^l, it is helpful in understanding the origin of the low frequency modes described 
by Eq.Q. This description also raises the question of negative energy modes, which 
can be supported in similar beam-plasma systems ^1]. Indeed, trapped electrons 
modes similar to the EAW described here can exist in a negative energy configuration 
|15|. although not in the regime considered here. 



EAW physics in LPI 



5 




Figure 2. The distribution (a) which exhibits local flattening at Vp, is 
mathematically identical to (b) , which comprises a background Maxwellian plus a 
smaller drifting electron distribution centred at Vp . The Doppler shifted frequency 
of plasma oscillations supported by the drifting electron distribution governs the 
EAW dispersion relation in the limit of small density. 



fedv - (8) 



3. Simulating an EAW 

A full treatment of the EAW requires a kinetic description of the plasma. This 
section outlines the fully kinetic Vlasov-Poisson model and the development of initial 
conditions for, and the simulation of, a travelling EAW. The model used is a 
one-dimensional Vlasov-Poisson system of electrons and immobile protons with no 
magnetic field in a numerical implementation which has been used previously 
to explore kinetic phenomena relevant to laser-plasma interactions This fully 
nonlinear self-consistent system is governed by the Vlasov equation for the electron 
distribution function /e, 

^ + ^?r--^§^-o (7) 

at ox me ov 
and Poisson's equation for the electric field 
dE__ e_ 
dx eo 

As an initial condition, let us consider an unperturbed distribution function fiattcncd 
at a phase velocity Vp, such that Vp = uj/k where oj and k can be chosen from the 
EAW branch of the dispersion relation. We specify 

/u = /o + /i (9) 

where /o represents a Maxwellian distribution and /i is given by Eq. © . The width of 
the flattened region, which relates to the number density of trapped electrons, is given 

1/2 

by At; and is proportional to the EAW amplitude (explicitly. Aw = (eEo/knie) ' ). In 
order to create a travelling wave we perturb the distribution function given in Eq.® 
by setting /e = /« + fp where 

fp{x,v) ^ 7 T—-sm(kx)dyfu (10) 

■me(uj — kv) 

Equation (|10|l represents a finite-amplitude generalisation of a linear perturbation 
of the Vlasov equation about /„. Note that Ec ljlUI) contains no singularities since 
dvfu\vp = and uj is chosen to be real. The Vlasov-Poisson system is initialised 
against a neutralising ion background in a periodic box with the distribution function 
fe = fu + fp and ui = O.Guipe, k = OAXJ-,^, Sn = O.lUe where 6n — Eq/L^- In this 
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Figure 3. (a) Surface plot of the electron distribution function for a large 
amplitude EAW (<5n = 0.1ne,u = 0.6ujpe,k = 0A\^^) at time t = lO^ujpc, 
simulated using the Vlasov-Possion code. (b) Corresponding contour plot. 
Contours for < 0.135 are drawn with dashed lines to highlight the trapped 
electron phase space structure. The v axes are given in units of VTe and the x 
axes in units of A/j. System boundaries are periodic. 
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Figure 4. Logarithmic electric field amplitude of EAW (5n = 0.1ne,c<j = 
0.6u!pe,k = 0.4/Aj3) against time. After an initial transient stage, the EAW 
persists as an electron plasma wave with frequency below the plasma frequency, 
which is undamped except for limited damping due to numerical diffusion. 



regime of {lu, k) we would not normally expect an electron plasma wave to propagate, 
undamped or otherwise. Figure 13 shows the trapped electron distribution of the 
EAW after a thousand inverse plasma angular frequencies. Figure 0] shows that the 
amplitude of the EAW is effectively constant: after an initial transient phase, only 
weak numerical damping remains. This numerical approach thus demonstrates how 
a non-Maxwellian distribution, specifically the flattening at the phase velocity of the 
wave introduced through Eqs.® and is necessary for the propagation of an EAW. 
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4. Relevance to the Langmuir Decay Instability 

A Langmuir wave can decay into a second Langmuir wave of lower wavenumber plus 
an ion acoustic wave (lAW). This process can occur repeatedly, forming a Langmuir 
cascade JS] . Can the EAW perform the role of lAW to produce a Langmuir cascade 
on electron timescales? The conventional Langmuir cascade proceeds for all k 
above a critical value determined by the point where the group velocity du/dk of 
the parent Langmuir wave (L) is equal to that of the lAW: 

dk ~ dk 3A, 

A similar analysis, for small /c, can be performed in the case where the lAW is replaced 
by an EAW. Approximating the Langmuir dispersion relation by 

w«Wpe(l + 3fc2A|,/2) (12) 

and the EAW dispersion relation by 

uj ~ LOpg^ (1.35fcAD) (13) 

gives a critical wavenumber kc ~ 0.45A^^, suggesting that LDI via the EAW might 
be a possibility. However, Eas. (|12(l and (|13|l are no longer valid for such a high 
critical wavenumber. The assumption of small k is therefore abandoned and the 
gradients calculated numerically to give Fig[Sl It follows from this full treatment, 
valid for all fc, that a process of Langmuir decay via the electron acoustic branch is 
not possible. However this does not rule out all forms of interplay between LDI and 
EAWs. The upper branch of the dispersion relation, essentially an undamped form 
of the conventional Langmuir mode, may replace one or both of the Langmuir waves 
in the LDI without affecting the critical wavenumber. This scenario is left for future 
work. 




5. Stimulated Electron Acoustic Scattering (SEAS) 

The collective scattering of incident laser light from an EAW can be simulated using 
a Vlasov-Maxwell model. Here the relativistic Vlasov equation for electrons, in the 
presence of transverse fields 

_|.^S)^ = (14) 
dt rUe dx nrie ^ ^ ^ dpx 

is solved against a stationary ion background, together with Maxwell's equations 

dEy _ ^dB, Jy 

dB, ^ dEy 

dt dx 

in one dimension. Transverse motion of particles is treated as fluid-like, hence 



(16) 



'-^=--Ey (17) 

dt rUe 

Jy = - erieVy (18) 

Poisson's equation (Eq.®) is solved as before, to give the longitudinal electric field. 

The initial conditions are chosen to prohibit conventional stimulated Raman 
scattering (SRS), so that density is above quarter critical density, and to satisfy 
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Figure 5. Difference between the gradients of the dispersion relations of the 
Langmuir and Electron Acoustic modes for a range of fc. In order for Langmuir 
decay to occur, this quantity must be greater than zero. The critical wavenumber 
kc is the point at which the gradients are exactly equal. This curve remains 
negative for all k, demonstrating that straightforward Langmuir decay via the 
EAW, rather than the lAW, is not possible. 



wavenumber and frequency matching conditions for SEAS. The system is periodic 
in X and undriven, with a TEM wave present throughout the system initially. Vlasov 
codes are inherently noiseless, so a low amplitude density perturbation is added to a 
Maxwellian velocity distribution to seed the growth of the EAW. The wavenumbers 
and frequencies of the incident wave El (fc^i = O.SlScjpeC^^, ujei — l-29a;pe), 
scattered wave E2 {kE2 = — 0.108wpeC~^, (jJe2 — l.OOSSwpe) and the EAW seed 
{kEAW — 0.923ujpeC''^ , ^EAW ~ 0.2842wpe) are chosen to satisfy the matching 
conditions for SEAS. The incident wave amplitude is Ei = 0.3e/{LOpeCme) and the 
system length ~ 30840^^^.^ . The simulation is conducted on a numerical grid with 
8,192 points in x and 512 in p. 

Figure IHl shows the electron distribution function at late time. The evolution 
of trapped electron structures, and resulting flattening of the distribution function, 
is visible, corresponding to an EAW. We thus have SRS-like scattering in a plasma 
whose density is greater than quarter critical: this scattering is from an EAW, an 
electron plasma wave with a frequency below the plasma frequency. 

This demonstrates the capability of a Vlasov based kinetic plasma model to 
simulate the scattering of incident laser light from non-Maxwellian, trapped particle 
distributions. The development of such instabilities in realistic parameter regimes, 
and in the presence of an externally driven TEM wave, is now considered. 
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Figure 6. (a) Surface plot of the electron distribution function at t = IQ'^ujp^ . 
Only a small section of the complete system is shown for clarity. Electron trapping 
and flattening of the distribution function can be seen: this is the EAW which has 
grown from a background density perturbation as a result of SEAS. Axes are given 
in relativistic units, c/i^ipe for space and m^c for momentum, (b) Corresponding 
contour plot. Contours below /e = 0.725 are drawn with dashed lines to highlight 
the trapped electron holes. The p axes are given in units of rrieC and the x axes 
in units of cwp^ . 



6. Kinetic effects in single hot-spot experiments 

Non-Maxwellian particle distributions, which require a Icinetic treatment of the 
plasma, can significantly affect the scattering of incident light, destroying the idealised 
picture of a three-wave parametric instability by allowing scattering from plasma 
oscillations omitted from conventional fluid treatments, such as the EAW. Section 
121 summarised the linear theory underpinning the EAW, section |21 demonstrated its 
existence in the non-linear regime by way of electrostatic Vlasov-Poisson simulations, 
and section [S] confirmed the possibility of stimulated scattering, resembling SRS, from 
an EAW in a plasma of greater than quarter critical density. We now utilise an 
expanded Vlasov-Maxwell code to investigate SEAS and related kinetic effects in a 
regime close to that achieved in single hot-spot experiments ^ 12] ■ This involves less 
than quarter critical densities (hence permitting SRS) and the presence of a continuous 
EM driver. 

6.1. Numerical Approach 

The Vlasov-Maxwell code, described previously, is expanded to allow for the presence 
of a continuous, sinusoidal, EM driver at a; = 0. This requires that the system no 
longer be treated as periodic, instead the boundaries are open. Any charge flowing 
past the system boundaries is assumed then to reside on a 'charged plate', external 
to the system. This external charge is included when calculating the electrostatic 
potential in order to avoid the creation of a DC field. The electrostatic potential (p is 
found using a tridiagonal matrix inversion, and the electrostatic field is then given by 
dxip- This replaces the Fourier method used to solve for in the periodic 
described in Ref . |TB] . 
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6.2. Initial conditions 

A system of normalised units is adopted in which time is normalised to units of uj~^ 
and velocities to units of c. Thus space is normalised to units of cw^e^ , electric fields 
to units of nif.cujpe/e and magnetic fields to units of meWpe/e- The laser intensity Jq, 
electron temperature Tg and density Ue achieved in single hot-spot experiments P| E] 
were, approximately: 

/o = 1.6 X IQi^'Wcm-^ (19) 
Te = 350eV (20) 
Tie = 1.2 X lO^Ocm"^ 0.03nc (21) 

These imply values for the simulation parameters (incident EM wave amplitude Ey^ 
and frequency loo, thermal velocity VTe and density n^) of 

Eyg = 0.33meCU!pe/e (22) 

ujq — 5.7775ujpe (23) 

VTe = 0.026c (24) 

He — 1 X ujp^eome/e^ — O.OSuc. (25) 

To minimise the charge loss from the system, a 'flat-top' density profile is used, where 
the density of both electrons and the neutralising ion background drops smoothly 
from no to zero over a distance « 40cuj~^ at the edges of the system. The simulation 
domain extends from a; = to x = 220cw~g^ , leaving a flat region at the centre of the 
simulation box approximately 140cw~/ in length, from p = — 0.75meC to p = 0.75meC. 
The simulation grid has 16,384 points in x and 1,024 points in p. The simulation 
runs to an end time of 1200u!p^. 



6.3. Results 

Figures [3 and |H1 display windowed Fourier transforms of the electrostatic field and of 
the back-propagating EM field, taken with a Hanning window of size « lbuj~^, at 
the centre of the system. These show the development of low frequency plasma waves 
after t = QQQLOp^. In the initial SRS burst, starting at t « AbQu)^^ the EM driver at 
Wo given by Ea. H23|) scatters from a Langmuir wave at loi = 1.06a;pe; k = 0.27/Xd, 
Vp = 3.93vTe, to produce reflected hght at a frequency uj2 = 4.72a;pe- This instability 
saturates via the trapping of electrons. Figure shows the electron distribution 
function during the late stages of the SRS burst, when electrons have been trapped 
and accelerated. A beam, similar to that observed in simulations of Raman forward 
scatterjSOj, forms in the electron distribution which is clearly visible in plots (FiglHt)). 
The trapping of electrons by the Langmuir waves driven through SRS evolves into 
a plateau in the electron distibution. This flattened region extends to low phase 
velocities, providing an environment in which low frequency plasma modes are able to 
grow and propagate. These low frequency modes are visible in the electrostatic field 
spectrum after the collapse of the initial SRS burst at t w 600a;~g^ , and correspond to 
two distinct electron acoustic waves (eawl and eaw2) at ujeawi — O.Tiujpe, k ~ 0.27 
Vp — 2.73vTe and later uJeaw2 — 0.57wpe, k = 0.28A2^^, Vp — 2.03vTe- 

The electron distribution at late times thus deviates significantly from a 
Maxwellian. The trapping of electrons in the initial SRS burst flattens the distribution 
around p = 0.1, allowing the development of low frequency plasma waves, the 
EAWs, whose trapped electrons further distort the distribution of particles. By 
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Figure 7. Windowed Fourier transform of the electrostatic field Ex at tlie centre 
of the system. An initial SRS burst at t 450/u)pe saturates via the trapping 
of electrons which distort the initially Maxwellian distribution and provide an 
environment in which waves below the plasma frequency can grow and propagate. 



The traces at ^ O.Sojpe and w ^ O.Gwpc, first appearing at t 600/a;pe, 
represent EAWs with phase velocities at Vp = 2.73fTe and 2.03vTe respectively. 
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Figure 8. Windowed Fourier transform of the backwards propagating EM field 
at the centre of the system. The spectrum shows the light scattered by Langmuir 
waves (SRS) and EAW waves (SEAS) identified in the electrostatic spectrum at 
the same point in space (see FiglTl. 
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Figure 9. (a) Surface plot of the electron distribution near the centre of the 
system at t = 500/ajpe. Electron trapping, visible here, is responsible for the 
saturation of the Raman instability and the creation of the electron beam in the 
spatially integrated distribution, (b) Spatially integrated electron distribution 
functions, for t = 0, SOOojpe and lOOOiJpe, normalised to the initial Maxwellian 
distribution. The trapping of electrons in the Langmuir wave driven by SRS 
temporarily creates a beam structure. The collapse of this structure is responsible 
in part for the formation of a broad plateau in momentum space at late times, 
which supports EAWs. 

the simulation's end, it has become clear that the plasma, and hence the modes 
which it supports, is not well described by linear or fluid approximations. Scattering 
observed in single hot-spot experiments was from EAWs with phase velocity v = lAvTe 
(fc = 0.29A^^, u! = 0.41a;pe)j with a backscattered wave amplitude aproximately three 
thousand times smaller than that from SRS. The amplitude of EAWs, and of the light 
scattered from them, observed in simulations is greater than observed experimentally. 
The simulations presented here also produce EAWs with higher phase velocities (i.e. 
Vp « 2.7vTe and Vp « 2.0vTe compared to Vp « IAvtc) than the scattered spectra 
from experiments indicate. These two deviations are closely related. As shown earlier, 
the dispersion relation for the EAW is dictated in part by the mode amplitude. As the 
EAW amplitude is increased, the dispersion relation shifts inwards, as shown in Fig^ 
resulting in a higher phase velocity at fixed wavenumber. Further work is required 
to quantify in greater depth this inconsistency between numerical and experimental 
results. The simulation runtime, t = 1200ix>~g"'^, is equivalent to less than three 
picoseconds - this serves to highlight how rapid the switch from the fluid to the kinetic 
regime may be, at the laser intensities considered here. As laser intensity increases, 
the kinetic effects discussed here will become more critical to the understanding of the 
associated laser-plasma interaction physics. 

7. Conclusions 

Experiments studying fundamental laser-plasma interactions in a sing le hotspot [BE] 
observed backscattered light from the interaction of the incident beam with two 
distinct plasma modes. First, there is scattering from waves having high phase velocity 
Vp « 4:.2vTe and a frequency above the plasma frequency Wpe, which was attributed 
to SRS: the three-wave parametric instability involving a Langmuir wave. Second, 
there is scattering from waves of considerably lower phase velocity Vp « lAvxe, whose 
frequency is below w^e, as low as OAlujpe: these low frequency modes were identified 
as the electron acoustic wave. The simulations reported here have attempted to model 
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the key physics of this scattering using a ID Vlasov-MaxweU approach. These have 
been successful in achieving scattering from both high and low frequency electron 
plasma waves, but have not been able to reproduce exactly the phase velocities of the 
EAWs and the relative amplitudes of the scattering events. 

The electron acoustic mode is a counter-intuitive phenomenon with a sparse 
literature: an electron plasma wave which propagates, free from Landau damping, 
at frequencies below the plasma frequency. This work has sought to clarify its 
characteristics, in terms of dispersion relations and the role of electron trapping, 
which also present an interesting application of plasma kinetic theory. Accurate 
representation and evolution of the complete phase space is of importance to SEAS, 
and is also vital to the saturation of the Raman scattering instability and the 
subsequent evolution of the system, as demonstrated here. Our simulations, in a 
regime close to those achieved in single hot-spot experiments, highlight the importance 
of kinetic effects, and the effect that the evolution of non-Maxwellian particle 
distributions may have on the scattering of incident light from an initially homogenous 
plasma. Even when the energy flow associated with SEAS is small, SEAS may 
have future applications as a diagnostic of the electron velocity distribution, given 
theoretical understanding and an appropriate modelling capability. 

Recent work_21_ .22 , 23 has identified the need for a deeper understanding of laser- 
plasma interactions, particularly in the regimes currently being approached by the 
next generation of lasers. The accurate noise-free representation and evolution of the 
particle distribution functions provided by a Vlasov code make it a valuable additional 
tool complementing both fluid and particle-in-cell descriptions. While a full 3D 
Vlasov treatment is beyond the limits of current computing power, ID and 2D Vlasov 
systems are tractable and can address many relevant problems. The present study has 
indicated some of the distinctive features of the EAW and SEAS physics that arise from 
the fact that the EAW can be considered to be primarily supported by an energetic 
particle population. Recent numerical work has highlighted the possibility that 
'trains' of electron holes, with low phase velocities, may be created by the action of a 
strong electrostatic driver at a frequency above the plasma frequency. Such structures 
could become involved in SEAS, and may be excited by the electrostatic daughter 
waves driven by the stimulated Raman and Brillouin instabilities. The interplay 
between the electrostatic mechanisms outlined in Ref.[55 and the electromagnetic 
scattering mechanism outlined here is a potentially interesting topic for future work. 
Finally, some of the conceptual links to the role of energetic particle populations in 
magnetically confined plasmas were noted in the introduction, and these too may 
repay further investigation. 
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